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£Sj Higher spin gravity in three dimensions has explicit black holes solutions, carry- 
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ing higher spin charge. We compute the free energy of a charged black hole from the 
holographic dual, a 2d CFT with extended conformal symmetry, and find exact agree- 
ment with the bulk thermodynamics. In the CFT, higher spin corrections to the free 
£T) energy can be calculated at high temperature from correlation functions of W-algebra 
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currents. 
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1 Introduction 

Higher spin theories of gravity in Anti de Sitter space, with massless gauge fields of spin s > 2 
in addition to the graviton, are dual to field theories with extended conformal invariance [JJ [2] 
EHUEIE]. On the gravity side, higher spin gauge symmetries combine with diffeomorphisms, 
mixing the metric with the higher spin fields. Therefore a generalized notion of geometry is 
required in these theories, providing a simple toy model for string theory in the low tension 
regime where higher spin string states become massless. In the CFT, conserved higher spin 
currents dual to the bulk gauge fields impose powerful constraints that, in several examples, 
render the theory free or exactly solvable. 

In AdS3/CFT 2 , the duality relates higher spin gravity to CFTs with W-symmetry [3 [8]. 
In particular, 3d Vasiliev gravity [HI QUI [HI [12] is conjectured to be dual to the family of 
interacting but exactly solvable Wjv minimal models [13] . In the large N limit, the detailed 
microscopies of this proposal are supported by symmetries [2] , the spectrum of primaries 
[15] . and certain correlation functions [16], [T71 [181 H3 EO]; at finite N some new ingredients 
may be required in the bulk (or perhaps can be found within the Vasiliev theory itself [2~T]). 

In this paper we focus on a universal aspect of the correspondence: black holes in AdS3 
carrying higher spin charge. These black holes encode the thermodynamics of the CFT at 
high temperature, which in two dimensions is determined by the chiral algebra. Thus black 



holes are insensitive to some of the microscopic details of the CFT, but capture certain 
features of the symmetries, as well as the asymptotic density of states. 

It is not obvious how to define a black hole in higher spin gravity. The usual definition 
- a spacetime singularity hidden behind a horizon --is difficult to apply because neither 
the Riemann tensor nor the causal structure of the metric are gauge invariant. In Euclidean 
signature, the problem is simpler, because a black hole is simply a smooth classical solution 
with torus boundary conditions. This definition has been used to construct explicit black 
hole solutions carrying higher spin charge [22] (see also (23J EH EH]). An interesting feature 
of these constructions is that turning on a chemical potential for a higher spin charge means 
deforming the CFT by an irrelevant operator, so these black holes describe an RG flow to a 
UV fixed point. The flow can be solved exactly in the bulk, resulting in a new AdS geometry 
in the UV that is realized by a different subset of the bulk gauge fields. 

There are many different theories of higher spin gravity in AdS3, specified by a Lie 
algebra. The original black holes of |22j are solutions of sl(N) gravity, which has gauge 
fields up to spin N. We will consider higher spin gravity based on the infinite-dimensional 
Lie algebra hs[A], since this is the bulk theory in the minimal model duality. (It is also in 
a sense universal, because setting X = N gives sl(N) gravity.) Black hole solutions in this 
theory were constructed in [26] . The mass, angular momentum, and charges of the black 
hole were computed and used to infer the free energy [26] 
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where a is the chemical potential for spin-3 charge, f is related to the mass and angular 
momentum (see below), and c = 3£/2G with £ the AdS radius and G the Newton's constant^] 
The free parameter A is a choice of theory, labelling the higher spin algebra lis [A]. Our aim 
is to reproduce this formula from CFT; the usual thermodynamic relation also gives the 



entropy, S — (1 — (3dp)\ogZ. In [26J, the formula (1.1) was already derived from CFT at 
the special values A = 0,1, where the chiral algebra has a simple free field realization, but 
for arbitrary A a more general approach is needed. 

Black hole entropy was first derived from microstate counting in string theory by Stro- 
minger and Vafa, using the detailed microscopic description of D-branes [27]. Later, this 
computation was streamlined and generalized by showing that the asymptotic density of D- 
brane microstates can be counted indirectly using Cardy's formula, for which the only input 
is the central charge of the CFT [28]. We will derive a similar universal formula for the 
thermodynamics of higher spin CFTs, and show that it exactly matches the gravity result 



We display only the holomorphic piece of the partition function; the right-moving sector gives a similar 
contribution, in gravity and in CFT. 



0). 

The general strategy is as follows. We start with the partition function of the CFT, 



Z CFT = Ti{e~ m+2maWo \ , (1.2) 



where Wq is the zero mode of the spin-3 current. The CFT has Woo [A] symmetry, an extension 
of the Virasoro algebra by an infinite number of higher spin currents. The partition function 
is expanded perturbatively in a; then a modular transformation is applied to each term, 
which makes it possible to compute the trace at high temperature. Thus the higher-spin 
corrections to the partition function are accounted for by certain correlation functions of 
W's on a torus. Unlike the bulk calculation, where only hs[A] is used explicitly, the CFT 
result depends in an essential way on the nonlinear terms in Woo [A] . 

The computation relies only on the chiral algebra of the CFT, and does not require (or 
allow for) explicit enumeration of the actual microstates contributing to the entropy. This 
means that although the minimal model CFTs appearing in the duality proposed in [T3] 
certainly reproduce the black hole entropy at high enough temperatures, it is not a fine- 
grained test of the specific details of that CFT; we will return to this point in the discussion 
section. 

In section [2j we explain in detail the strategy of the computation, and state our main 
technical tool, a recursion relation for the evaluation of torus correlation functions. In section 
pi we apply the recursion relation and compare to the gravity answer (1.1). In section Hi 



we conclude with a discussion of the regime of validity of the formula and applicability 
to minimal model holography. Various technical details are relegated to the appendices: 
appendix [S] defines the Weierstrass functions that appear in torus amplitudes; in appendix 
[Bj we define some notation and derive the recursion relation for torus amplitudes; details of 
the six-point calculation are given in appendix [Cj and the commutation relations of Woo [A] 



appear in appendix D 



2 Higher Spin Corrections to CFT Thermodynamics 

Let us begin by explaining the aim of our calculation in more detail. The entropy calculation 
of [2B] amounts, from the dual CFT point of view, to the evaluation of the partition function 

Z C FT(?,a) = Tr(>>-£ y w °) , q = e 2 ™ , y = e 2 ™ , (2.1) 

where the trace is taken over the entire spectrum. Here Wq is the zero mode of the spin 3 
generator of Woo [A], with chemical potential a, while f is the complex structure of the torus, 



related to the black hole temperature Tjj and (imaginary) angular potential Qh by 

-tt H ). (2.2) 
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The bulk partition function (1.1) was obtained in [26] from the thermodynamics of the 
black hole solution in AdS3 higher spin gravity. It is an expansion in powers of the chemical 
potential a, so it should be compared to the CFT expansion 

Z CFT (r,a) = Tr(>-M) + ^^ Tr(> ) 2 $*•-&) + ^^ Tr(> ) 4 g^) + ... 

(2-3) 
in the high temperature regime, i.e. for t — >■ 0, and to leading order in the central charge 
c (since the gravity calculation is only reliable for large c). At high temperatures, the in- 
dependence of each term in the expansion is fixed by conformal invariance, which requires^] 

logZ C FT(r,a)wi/(^) (2.4) 

for some function /. 

As is familiar from entropy calculations |2S|, the standard method to obtain the par- 
tition function from a dual conformal field theory point of view is to do the S-modular 
transformation 

r = -I , q = e 2wiT . (2.5) 

T 

In the high temperature regime, i.e. for f — > 0, we have q — > 0. The answer for the trace is 
then dominated by the contribution from the vacuum state. 



This argument can be directly applied to the first term in the expansion (2.3), 



Tr(>>-s) =J2SsrTr r (q L °-^ ~ (j2 S °°) <T^ + ' ' ' , (2-6) 

s,r \ s / 

where the sum runs over all primaries labelled by r, s (with r = the vacuum representation), 
S sr is the modular S- matrix (not to be confused with the black hole entropy), and the dots 
indicate terms exponentially suppressed at high temperature. The leading behaviour of the 



Reintroducing the size of the system R, the dimensionful temperature T ~ i/{fR) and a chemical 
potential a ~ aR 2 of dimension —2, as appropriate for a spin-3 charge, the only dimensionless quantity 
is olJt 2 . Since exp(2maWo) defines, at least formally a unitary operator, the asymptotic behaviour of 
Zcft{t,ch) must be the same as for the usual trace, i.e. proportional to t _1 , and thus the structure must 



be as in (2.4 1. Alternatively, we may argue that \ogZ can be at most linear in the system size R since the 



free energy per unit size should be finite. 



logarithm is then 



logTr^-^J =__r + ... , (2.7) 



and this reproduces precisely the a-independent term in (1.1), using the relation r = —\ 



This is equivalent to the Cardy formula for the entropy. In the following we want to repeat 



this analysis for the other terms in the expansion (2.3). The main technical problem we 
need to understand is how traces with zero mode insertions behave under the modular S- 
transformation. 

Consider first the simpler case where the zero mode is that of a dimension-1 current, and 
the modular transformation of the full exponential is actually known to be described by a 
Jacobi form. (This results in a simple Cardy-like formula for the entropy of a black hole 
carrying U(l) charge.) Let us assume that we have a rational CFT with a U(l) current J n , 
and define the trace in the representation labelled by r 

(j>r{r,z) = Ti r (y Jo q L °-^) , y = e 2niz , q = e 2niT . (2.8) 

Then it follows from usual CFT argument qj that 

, f(LT + b z \ v-^ r , cz 2 -] , , i a b\ , . 
*r (^. ^r d ) = E -p{ 2 ««^} «™ ( e d ) Mr, ,) , (2.9) 

where the sum runs over all irreducible representations of the chiral algebra, and M rs defines a 
representation of the modular group. Furthermore, the level k is defined via the commutator 

[J m , Jn] = 2km5 m _ n . (2.10) 

Expanding in powers of z, one then finds for example that under the ^-modular transforma- 
tion 

Tr r (Vo<f°-^) =^^[r 2 Tr s (j J g L °-^) +^rTr s (g L °-^)] . (2.11) 

s 

In the high temperature q — > limit, it is the second term evaluated for s = that is 
dominant. The zero modes inside the first trace annihilate the vacuum so it is exponentially 
suppressed. 



Unfortunately, the analogue of (2.9) for zero modes of general spin is not known. There 



have been some ideas that the relevant transformation property will be of the same structure 

3 In particular, this argument relies on the spectral flow automorphism of the U(l) current algebra, see 
for example [29] . 



as that described in [30], see [31] . but no general formula is knownn As we shall explain in the 
following, we can derive the transformation properties of these traces from first principles. 

2.1 The general strategy 

In order to explain the general strategy let us introduce a little bit of notation. For each 
representation of the chiral algebra (labelled by r), we can define a torus amplitude by 

F r ((a\ Zl ),..., (a n , z n ); r) = z 1 ? 1 ■ ■ ■ z h n n Tr r (v{a\ z x ) ■ ■ ■ V(a n , z n ) q Lo ~^ , (2.12) 

where hj are the conformal dimensions of the chiral fields a J , i.e. L^ = h^a? with hj G N. 
(In our case, the chiral fields will be the higher spin fields of the Woo [A] algebra.) These 
functions are periodic in each variable Zj under the transformations 

Zj i—)- e 2m Zj , and Zj i— > qzj , (2-13) 

and therefore deserve to be called 'torus amplitudes'. The periodicity under Zj H- e 2m Zj is 
obvious; to see the periodicity under Zj H- qZj we first note that the functions are independent 
of the order in which the vertex operators are inserted (since they are local fields). Using 
the transformation property of the fields under scaling 

V(a j ,q Zj )q L °-^ = q~ h i q L °~^ V(a j , Zj ) , (2.14) 



as well as the cyclicity of the trace it follows that the functions in (2.12) are also invariant 
under Zj h-> qzj. 

The amplitudes have a simple transformation property under the modular group provided 
that all a- 7 are Virasoro primaries so that L n o? = for n > 0. If this is the case then under 
a modular transformation 



F r ((a^^ 1 ),...,(a^^;f3) = (cr + rf)^^^M rs F s ((a^^^ +d ),...,(a^^ +d )^), 

(2.15) 
where M rs is again a matrix representation of the modular group, i.e. it depends on the 
matrix (" ^), but not on the a 1 . Note that this transformation rule makes sense: the left- 
hand-side is invariant under 

Zj h-> e 27Ti Zj , and Zj H- qzj (with q = e 2 "^) , (2.16) 



4 We thank Jan Manschot for bringing these papers to our attention. 



under which z ( j T+ transforms as 



z- ^ e z j 



and 



„CT+d l __, 2iri(aT+b) „cr+d 

z j \-> e z j 



(2.17) 



respectively. These are both periodicities of the right-hand-side, see (2.13) 



We are interested in the modular transformation properties of the traces with the insertion 
of zero modes. Expanding the vertex operators in modes as 



V(a,z) = J2 



&m. Z 



-m—h 



(2.18) 
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where h is the conformal dimension of a, the zero modes are obtained from the torus ampli- 
tude via the integrals § — , i.e. 



Ti r [a 1 ---a^q L °-^ 



(2ttz) 



dz\ 



dz n 



F r ((a 1 ,2;i),...,(a n ,z n );r) 



(2.19) 



Since we know the modular transformation properties of the torus amplitudes, this can now 
be used to deduce that of the traces with zero modes. Concentrating on the ^-modular 



transformation that is of relevance for us, we thus obtain from (2.19) and (2.15) 



(2m) n Tr r (a* ■ ■ ■ a n q L °~^ = ]T S rs r~ n+ ^ h > J 



1 Av 



dz\ 

Z\ 



dZr, 



Fs({a\z x ), 



,(a n ,z n );r) 



1 ^n 



(2.20) 
Here the integration contours are the result of the transformation z — > z T = z, which swaps 
the two cycles of the torus under the ^-transformation, see figure [TJ The contours are 
displaced to avoid short-distance singularities, so the contour for Zj runs from Ej to Sjq for 
some arbitrary phase Ej\ these phases are implicit throughout. From now on we work on the 
z-plane and drop the tildes. 



In the high temperature limit f — > the dominant contribution to (2.20) will come from 



the vacuum s = 0. Thus our calculation is reduced to doing the multiple integrals of the 
torus amplitudes in ( 2.20| ) with s = 0. In order to evaluate the torus amplitudes we can use 
the recursion relation of [32J to rewrite (from now on we drop the index s = from the torus 
amplitudes and the traces) 



F [(a 1 , zx), (a 2 , z 2 ), . . . , (a™, z n )\ rj 

= 4 2 ■ ■ ■ z h n " Tr (a} V{a\ z 2 ) ■ ■ ■ V(a n , z n ) q L °~^ 



(2.21) 
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Figure 1: Contours of integration for the insertion of two zero modes, before 
and after the S modular transformation. The fundamental region of the torus lies 
between the dashed lines, \q\ < \z\ < 1. 



where V m+ i(x,q) are the Weierstrass functions defined in appendix [Al and the bracketed 



modes, defined in appendix B.l have the form 



a[m] = (27ri)- m - 1 ^c(/i a 



,i,m)a- ha+1+i 



(2.22) 



Note that both sides are periodic under Zi i— > qZi, this is manifest for all terms except for 
the term in the second line, and the m = terms in the double sum of the last line, given 



that V\ transforms as (see eq. (A.4)) 



p[{ ^ q ) =Vl W q] 2 " 



(2.23) 



It is not difficult to check that the non-periodic contribution of (2.23) cancels precisely the 



non-trivial term that appears from the second line, where we pick up a commutator term 
from moving V(a l ,Zi) past a^. Indeed, this commutator is precisely equal to 



[a], V{z Lo a\ z)} = (2m) V(z L ° a 1 [0]a\ z) 



(2.24) 



so that 



F (. . . ,[a ,V(a\z)}, . . .) = (2m) F (. . . ,V(a[0]a\z), . . .) . 



(2.25) 



Unfortunately, this simple recursion relation is not quite sufficient (for n > 2) since 
the term in the second line involves an explicit zero mode, and therefore cannot be recursed 
further in this manner. We have therefore generalised the analysis of [32] to find also recursion 
relations for correlators involving zero modes. Defining 

F (bi; (a 1 ,z 1 ),..., (a n , z n ); r) = # ■ ■ ■ zfr Tr (V V(a\ z±) . . . V(a n , z n ) g Lo ^) (2.26) 



we have 
F^ia^zJ, . . . ,(a n ,z n );r) = F (blal-,(a 2 ,z 2 ), . . . ,(a n ,z n );T) (2.27) 

+ E E E ( £ ) 9Ui (-) F K\ ifl\ Si), • • • , {d^[m]a\ Zj ), . . . , (a«, z n ); r) , 



i=0 j=2 meN 



where 



and 



<4 +1 (x, g) = (2tt»)^ ^l d l T V m+1 ^x, q) (m > i) (2.28) 



?/;.! 



d (i) = (-i)*^])^ 1 . (2.29) 



This identity is proven in appendix |B|p| Note that the sum over m is actually finite because 
modes with high enough mode number will annihilate a? . 



With the help of this recursion relation we can then evaluate the amplitudes in (2.20) 
explicitly (and recursively). The recursion formula turns the amplitude into a sum of prod- 
ucts of Weierstrass functions; then the contour integrals are straightforward because V(z)/z 
is a total derivative. This computation will be explained in detail for the (Wo) 2 , (Wo) 4 , 
and (IVo) 6 cases in the following section. Fortunately, we do not need to evaluate the exact 
answer, but only the leading order contribution as c — > oo. 



3 Applying the Recursion Relation 

Exponentiating the gravity result (1.1) leads to 

-c/24 



1 n . / 1 2 4 50 A 2 - 7 4 8 200 5A 4 - 85A 2 + 377 , 12 \ 

(3.1) 
where we have written only the order c part of the expansion. Higher powers in c are given 



by disconnected diagrams, so we will compare (3.1) to the connected contribution to the 
CFT partition function. To compute connected correlators, we simply discard terms where 
some subset of the W's produce a central term at an intermediate step in the calculation. In 
this section we will use the recursion relation to compute correlation functions of W currents 



and thereby reproduce (3.1) in CFT, term by term. 



5 The function g l m+1 (x, q) is extended to the range m < i in appendix lAl This is necessary for the general 
recursion relation if zero modes of the chiral algebra do not commute, Dut will not appear in any of our 
computations. 



3.1 The 2-point case 



The first higher-spin correction to the partition function is independent of A, and was derived 
from CFT in [26] by taking advantage of the free-field realization of Woo [A] at A = 0, 1. We 
will rederive it by the strategy outlined above, in order to illustrate the method. 



The leading correction in the expansion (2.3) is, using (2.20) 
z {2) = (2maf 



„ Tr((^ )V°--) -^f- r d ^F((W,z 1 ),(W,z 2 );r) . (3.2) 
The symbol '~' means we have dropped terms that do not contribute to \ogZ at leading 



order. Applying the recursion relation (2.21) 



F((W,z 1 ),(W,z 2 );t) = zlTi(WoW(z 2 )q L °-^) +V m+1 (^j F((W[m]W,z 2 y,r) . (3.3) 

Here, as in the following, we shall always imply a sum over the Weierstrass index m > 0. In 
the first term, only the zero mode of W(z 2 ) survives in the trace, but Tr ((Wo) 2 q L °~ c ^ 24 ) has 
no contribution from the vacuum state and hence is exponentially suppressed. In the second 
term, the vacuum contribution to the trace dominates, so 



* J l z l J l z 2 \ Z l. 



(3.4) 



where (•) is the vacuum expectation value (on the sphere). The term with m = vanishes by 



(2.25). According to (A. 2), integrating a Weierstrass function reduces its index by 1. Only 



V\ is not periodic, so this implies 



q dz 2 _ I ' Z\ 

— v 2 (- 

i z 2 \z 2 



(2tu) 2 , 



1 z 2 



2 i-i I Z i 

Z 2 



(m>l) 



Therefore 



Z (2) w ± q -c/U(2<iri)* a 2 r i (W[l]W- 3 ) 



The definition of the bracketed mode (B.5[) is 



W[l] = (2iri)~ 2 W-i + -W + -W l W 2 + —W 3 + 
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30 



(3.5) 



(3.6) 



(3.7) 



Since we are only interested in the vacuum expectation value, the only surviving term is 



then (W[1]W_3) = (2m) -^{WsWs), and plugging in the commutators of Woo [A] from 



10 



appendix D the final answer is 

Z (2) =2mc ^_ g -/24 ) (3i 



aV 



in precise agreement with the leading correction to the gravitational partition function (3.1 ). 



3.2 The 4-point function 

Now we proceed to the a 4 correction to the partition function, 

Z {i) = ^^Tr((^ ) 4 ^-^) ^^ j F((W,z l ),(W,z 2 ),(W,z 3 ),(W,z A );r) , (3.9) 
where throughout this subsection, 



dz\ f q dz^ 

i z i J\ z k 



(3.10) 



Applying the recursion relation (2.21), we get 



F((W,z 1 ),... J (W J z 4 y,r)=JF(W ;(W,z 2 ),(W J z 3 ),(W,z 4 y,r) (3.11) 

+ 3jv i+1 (^)F((W,z 2 ),(W,z 3 ),(W[£}W,z,);T) , 

where sums over Weierstrass indices are always implied. Three terms on the right have been 
combined by relabelling integration variables; strictly speaking, this is not allowed, because 
the integration contours are displaced by phases Sj as in figure [ll but we will account for 
this subtlety below. Note that we cannot drop the first term on the right, despite the zero 
mode which naively implies that the leading 0(q~ c / 2A ) contribution vanishes. The reason is 
that q also appears in the limits of integration, and after integration we will see that this 
term contributes at leading order. 



11 



Applying the recursion relation twice more to (3.11) gives 
JF((W,z 1 ),...,(W,z i y i r) = 



SjV i+1 I j\ V m+l I Ij V k+ i (j) (W[k]W[m]W[£]W. 3 ) 

+ 3JV W (A V m+l (A V k+l [A ((W[m]W)[k)W[£)W- 3 ) 

+ 5 -^jVe +1 (A d T V m (A (dV[m]W[t]W^) (3.12) 

+ ^ fdrVt (A -Pm+1 (A (W[m]dU[e\W„ 3 ) 

+ (2«) 2 



Z2j \Z 3 

dlVl-i (A (d {2) m-z) 



where (•) is the connected part of the vacuum expectation value and we have defined the 
states 

rfW = -W[0]W , and d (2) = W[0]W[0]W . (3.13) 

Each of these integrals can be computed easily using properties of the Weierstrass functions, 
see appendix [A] for details and examples. Some of the integrals diverge if we ignore the phases 
£j that separate the integration contours. For example, making the contour separations 
explicit, the last line includes 

Hr^ffe^W^f?!), (3.14) 

Je 3 z 3 J £4 z A \z 3 J \e 3 J 

which diverges as £4 — > e 3 . However, since the original expression is independent of the 
contours, we can account for all of the contour dependence by simply dropping divergent 
terms of the form V m (l). These terms must cancel in the sum, and the pole structure (as 
functions of the Ej) requires that they cancel only with other Weierstrass functions of the 
same order, so they do not leave any finite remainderjj 

The same argument implies that we can drop terms proportional to r 2 or higher powers 
of t. As noted in section [2j the powers of r and a in the partition function (1.1 ) are dictated 



by conformal invariance. Only the r 1 term in (3.12) has the correct power; other powers 
can only appear when multiplied by appropriate ^-dependent functions (or suppressed by 
powers of q), and such terms cannot contribute. 

As a consequence, only terms with bracketed indices summing to 3 need to be computed 

6 As a check, we have confirmed by explicit computation that the contour-dependent terms in the 4-point 
recursion cancel in the final answer. This requires keeping track of contours at each step in the recursion 
and computing many additional terms. 

12 



in gig). That is, 



y F((W, Zl ), ...,(W, z 4 ); r) = 3/!^ + A 2 ) + 3/ 2 A 6 + ^/ 3 A 3 + ^ 4 + hA 5 , (3.15) 
where the integrals are 



h 

h 



(27Ti) 6 

1 



*'>£K 



P2 ,|^2 



(27T«) 5 

and the algebra factors are 

A x = (27ri) 6 (W[l}W[l]W[l]W_ 3 ) 
A 3 = (2mf(d^[2]W[l]W. 3 ) 
A 5 = (2m) 6 (W[l}d^[2}W. 3 ) 



h 

h 



(2mf 

1 
{2mf 



*'>©*£ 



*<th*e 



(3.16) 



A 2 = (27T2) 6 ((^[1]^)[1]^[1]^_ 3 ) 

A 4 = (27rz) 6 (rf (1) [3]^[0]iy_ 3 ) (3.17) 

A 6 = (2m) 6 {(W[0]W)[3]W[0]W. 3 } = -A 4 . 



Here the identity between A 4 and Aq follows from the definition of d^\ To compute the 



vevs in (3.17), it suffices to ignore nonlinear terms in the W-algebra, as they contribute at 



subleading order in 1/c (though this will no longer be the case in the six-point function). 
Furthermore, we discard disconnected contributions, as described above. Nested bracket 
modes like (W[1]W)[1] can be computed directly by evaluating the state W^[l]jy_3f2, com- 
puting the corresponding [l]-mode, etc. A perhaps simpler approach, especially for higher 
n-point functions where nonlinear composite states contribute, is to use the Jacobi identity 
(see 



{aim 



]*)[*] = £ 



i=0 



• / m 



(-l)M . )a[m-i]b[n + i\- (-1) 



rn+i 



m 



b[m + n — i]a[i] 



(3.18) 



to remove the nested modes. Applying this identity, plugging in the definition of the bracket 



modes (B.5), and using the Woo [A] commutators, we find 
A x = c 



Ni 26iV 4 



9 135 



Ao = c 



4JVf , 2iV 4 



+ 



27 15 



(3.19) 



and 



A 3 = A 2 , A 4 = -A, = -A 2 , A 5 = A ± . (3.20) 

The Weierstrass integrals A, 2,3,4 are computed in appendix IA1 Plugging in the answers it 



13 



then finally follows from (3.9) and (3.15) that 



Z& « _^/24 27rzc A(5 A r2_7Ar 4 ) a 4 r 9 (3.21) 

« -g-^vnc^l^V (3.22) 



in precise agreement with the gravity answer (3.1) 



3.3 The 6-point function 

A new ingredient appears in the a 6 correction to the partition function: nonlinear terms in 
the algebra contribute at leading order. Schematically, the spin-3 and spin-4 commutators 
of Woo [A] have the form 

[W,W] ~ U + L + -L 2 + c 

c (3.23) 

[U,W] ~ X + W+ -LW , 

c 



where U is the spin-4 current and X is the spin-5 current (see appendix D for details). 
Note that nonlinear terms always appear with powers of 1/c, which is what allowed us to 
drop them in the 2- and 4-point calculation, but starting from the 6-point case they can 
nonetheless contribute. For example, starting with the the 6-point function we can contract 
two currents to make L and two currents to make U, 

WWWWWW -» WWW . (3.24) 

Now contracting UW ~ -LW + X + W gives the nonlinear term 
i i c 

-LLWW . (3.25) 

c 

Since LL and WW both have central terms, this term is of order c, and therefore contributes 
to log Z at leading order. 

Otherwise, the computation of the 6-point function proceeds as above. Starting from 



r(6) =- ( 27ria ) 

' ~ 6! 



Tr((^o)V "^)~^ r yi 71 ((^,^i),...,(iy,^ 6 );r) , (3.26) 



we apply the recursion relation to reduce the six-point correlator to a sum of vevs, perform 
the Weierstrass integrals, discard contour-dependent terms (i.e. V m (l) and terms not linear 
in r), and compute the sum. Some additional details are provided in appendix O The final 
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result is 



y V 648 972 ° 12150iV 3 2160 y V ; 

„ /94 200 5A 4 -85A 2 + 377 fi „ , 

« g- c/24 2vrzc— ^-^ «V 3 (3.28) 



in precise agreement with the black hole value (3.1). 



4 Discussion 

As we have shown in this paper, we can reproduce the higher spin corrections to the black 
hole entropy from integrals of correlation functions of W-currents on the torus. This was to 
be expected on general grounds, but it is satisfying to see that things work out correctly. In 
particular, the CFT calculation depends on the full non-linear Woo [A] structure and hence 
our result provides a detailed check that Woo [A] is indeed the correct symmetry algebra of the 
dual CFT. We should also stress that, on the face of it, our computation is entirely different 
from the bulk calculation performed in [26]. In the bulk, the entropy is determined by the 
solution of a zero-holonomy condition for hs[A] Chern-Simons gauge fields on a torus. The 
zero-holonomy condition ensures smoothness of the Euclidean horizon, which determines 
the thermodynamics of the black holejj Somehow, the holonomy condition encodes the 
full nonlinear structure of Woo [A] in a nice geometric way; it would be very interesting to 
understand how the Chern-Simons calculation emerges directly from CFT or vice versa. 



Let us close by commenting on the regime of validity for the free energy formula (1.1) 
and its interpretation in higher spin holography. In any unitary CFT with Woo [A] symmetry, 
assuming the theory has a gap, the formula applies asymptotically as T — > oo. However, 
this fact alone is not enough to conclude that a CFT with Woo [A] symmetry reproduces 
the thermodynamics of higher spin black holes. Indeed any CFT obeys the Cardy formula 
asymptotically, but not every CFT describes black holes in AdS3. One key feature that 
distinguishes holographic CFTs from the rest is that the Cardy formula applies not only 
asymptotically, but at temperatures down to the Hawking-Page transition in the bulk. These 
are very special CFTs, like the D1-D5 system dual to string theory in AdS3, with a large gap 
in dimensions separating excited string states from the low-energy fields of Einstein gravity, 
see e.g. 



For the case at hand, the dual CFT must have free energy (1.1) whenever the black hole 



dominates the bulk thermodynamics. However, it is not known whether Vasiliev gravity in 



7 Using this method, the authors of [53] were able to compute also the a s term in ( 1.1 ), whereas to do so 
in CFT by our method would be very difficult. 
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three dimensions has a Hawking-Page transition, or whether the black hole dominates the 
bulk thermodynamics anywhere besides T — > oo. If there is indeed a phase transition above 
which the black hole dominates, then the dual CFT should have a gap large enough so that 
(1.1 ) applies above the transition temperature. The microscopic CFT proposed in [13], a WV 



minimal model, has a large number of light states with dimension A < 1, so it presumably 



obeys (1.1) only at asymptotically high temperatures. This may be mirrored by the fact 
that the Vasiliev gravity theory has other, black-hole-like solutions [2T] which contribute to 
the bulk thermodynamics. 
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A Weierstrass functions 

The Weierstrass functions are defined by the power series 

(2m) k >r^ fn k ~ 1 x n 



^H(H)T£(j^), (a.! 



for A; > 1, which converge for |g| < |x| < 1. (This differs slightly from the standard 

Weierstrass function p k , see [32] for the relation.) They satisfy the important recursion 

relation 

d k 

We can rewrite V\ (x, q) as 

2iri y^ ( %n< f 



n^O 



?,(*,,) = ^ - 2« + 2™ JJ [f^- , (A.3) 



thus defining a meromorphic continuation of V\ (and hence, because of (A. 2), for all Vk) to 
|g| < \x\ < \q\~ x - A straightforward calculation then shows the identity 

Vi(qx, q) = Vi(x, q) + 2m , V k (qx, q) = V k (x, q) (k > 1) . (A.4) 



16 



In the recursion formula (2.27), we also need the closely related functions 

(2ni) k 



9k(x,q) 



(k-l)\ 



J2n k - i ~ 1 x n d t T 



n^O 



T \i-q n 



Note the relations 



9k{x,q) = V k (x,q) 



and 



gl(x, q) = (2m) i {k n % - -^ & T P^(x, q) (k - i > 1) . 



(k-l)\ 



(A.5) 

(A.6) 

(A.7) 



A.l Weierstrass integrals 



Products of "Pfc and g k can be integrated using (A.2) and (A.7). Dropping from now on the 



argument q, we have for example from (|A.2|) and (|A.4|) 



q dz, 



V 2 

1 Z 2 \Z 2 



(2ni) 2 , 



q dz. 



1 z 2 



2 v ' Zl 

I m+1 



'-2 



(m>l) 



(A. 



The integral of V\ can be computed by relating V\ to the Weierstrass (^-function (denoted 
pi in [32]), which integrates to the log of the Weierstrass a-function, giving 



q dz> 



2 v 1 ' Zl 



(2iri) (in — 2nir + log z\) . 



(A.9) 



'1 z 2 \^2y 

Other useful identities can be derived by integrating by parts, for example 

9 — logzt V m+1 (-) = (2m) 2 - V m (-) - (2m) 3 -S m , 2 (m > 2) . (A. 10) 

i z x \z 2 J m \z 2 J m 

We will also need to integrate g l k (x,q). In all of our computations, k > i + 1, so we can 
replace g\ (up to a constant factor) by d^Vk+i-i- Integrals of d l T V m can be computed by 
differentiating the identities above. For example, 



q dz 2 lz\ 



d T / —V m (-\ =0 (m>\) 
h Z 2 \z 2 J 

implies, acting with the derivative on the limits of integration and on the integrand, 



(A.ll) 



O t I m 

1 z 2 \z 2 



-(2ni)V m ( Zl ) (m>\). 



(A.12) 
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Differentiating again, 

T — d 2 T V m (-) = -(2m) 2 d T V m (zi) - (2m) mV m +i(zi) (m > 1) . (A. 13) 

Jl z 2 \Z2/ 

Combining these techniques allows one to compute all of the integrals encountered in the 



partition function calculations. For the four-point function, the relevant integrals from (3.12) 



are listed in (3.16). I\ is computed by integrating first over zi,z 2 ,z 3 , then over z 4 , which 
immediately gives 

h = 2mr . (A. 14) 

To compute I 2 , we first integrate over Zi,z 2 , and then integrate by parts to get 

I 2 = \ (2mr) - \r 2 V 2 (l) ->■ \ (2mr) , (A.15) 



where in the last step we have used that the divergent term P 2 (l) cancels if one carefully 



accounts for contour separations (see section 3.2). 



For ^3, we integrate first over z\ then apply (A. 12), resulting in 



I 3 = -2irir . (A. 16) 

Similarly 1^ is determined by first integrating over zi,z$, and then integration by parts, 
leading to 

h = \ (2mr) - \r 2 V 2 {\) -> \ (2mr) . (A.17) 

All of the integrals encountered in the 6-point function can be computed similarly. 



B Torus recursion 

B.l Bracketed modes 

Let us define the expansion coefficients c(h,j,m) via 

(log(l + z)) s (1 + z) h - 1 = J2 <h,j, s) J . (B.l) 

Then the bracketed modes can be defined using the 'Mathematician' convention for modes 
as 

oo 

a[s] = (2ni)- s - 1 Y^ c(h a , i, s) a(i) . (B.2) 



Here we have assumed that a has definite Lq eigenvalue, L^a = h a a, and the modes a(i) are 
related to the modes Oj via 

a(n) = a n+1 _ ha (B.3) 



since we have 



V(a,z) = ^a{n)z- n - 1 = J2 a nZ~ n ~ ha = J2 



dn+l-h a Z 



-n-1 



(B.4) 



Thus in terms of the usual physicists' modes we can rewrite (B.2) as 



a[s] = {2m) s l ^ c(h a , j + h a - 1, s) aj . 

j=S+l—h a 

For the Virasoro modes a different convention is usually used, namely 
L [n] = (2m)- n Y, c(2,j,n+ !)%_! - (2™) 2 ^ n ,_ 



j>n+l 



(B.5) 



(B.6) 



The relation to the above definition is Lr n i = {2ni) 2 {uj — £l)[n+ 1], where u is the field with 
mode expansion 

V(u, z) = J2 u{n)z~ n - 1 = J2 L n z~ n ~ 2 ■ (B.7) 



Note that (u> — ^1) has definite Lp] eigenvalue since 

L[Q\ = Lq + 2-^1 — g-^2 + • • • 



and hence 



L [0] (L_ 2 fi - £fi) = 2L_ 2 ft - |L 2 L_ 2 fi = 2(L_ 2 ft - £fi) . 



(B.8) 



(B.9) 



In these conventions we then have the identity 

(Lr_i]o)[n] = -na[w — 1] , 

which becomes for states of definite Lpi eigenvalue the identity (Lr_iio)r n i 
see page 4 of [33] . 

B.2 Derivation of recursion relations 



(B.10) 

■([h a }+n)a [n] , 



In this appendix we prove the recursion relation (2.27). The case without zero modes is 



derived in [32] and reviewed in [31], and a similar strategy is followed here. We begin by 
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expanding the first vertex operator in modes, 

F(b e ; (a 1 , z x \ . • • , (a n , z n ); r) = F (b e a ; (a 2 , z 2 ), . . . , (a n , z n ); r) (B.ll) 

We will now cycle a\ around the trace to evaluate the second line. For this we need the 
commutator 

[alVia^Zj)} = JT ( Hl ~^ + k ) V(al_ hi+1 ai,z,)z^ 1+k ~ m , (B.12) 

m=0 ^ ' 

and the relation 

a \ q Lo ~^ = q k q L o-M a \ . (B.13) 

When we move a\ past the zero modes b e we will pick up the commutator [o^,6q]. Let us 
define the modes of the single commutator to be d k — — ^[^o, ofy, and more generally, 

4 0) = 4 , 4° = -^-^ d t l) \ (i > 0) - (B.14) 

This definition can be recast in terms of states 

d {i) = (-I) 1 (b^Ya 1 , (B.15) 



as follows from (2.24). Next we prove the identity 



^ k F(b e al;(a 2 ,z 2 ),...,(a n ,z n y,r) = (B.16) 



A^m^/v^-/ ! 



^W-UJSw 



£ C * + fc ) F ^ ( ° 2 ' 22) " " ' {d *- h+iaj > Zj) > ■ ■ • ' (an ' Zn); r ) 



00 

X 

t=0 



where in the last line, h is the dimension of S l >. (Typically $- % > is not homogeneous under 
Lq, in which case this expression is defined by expanding in Lq eigenstates.) This is proved 
by induction on £. Moving a\ around the trace on the l.h.s., we pick up a commutator with 
each vertex operator and a commutator with b e . When it returns to its original position, the 
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trace is multiplied by q k from (B.13). Solving for the original expression, 



z 1 k F(b e a 1 k -(a 2 ,z 2 ),...,(a n ,z n );r) 



1 _ Q k 1^ [ z j 1^ 



" . \ v^ fh\ — 1 + k 



q*. a — i \ z 

' 3=2 V 1/ t=0 



(B.17) 
F (6j; (a 2 , 22), . . . , (a]_ hl+1 a 3 , Zj), ..., (a n , z n ); r) 



-fe 9 



' x l-g< 



: F([b i ,al];{a 2 ,z 2 ),...,{a n ,z n );T) . 



The second line is the i = term in (B.16), which proves the identity for £ = 0. To evaluate 



the third line, we plug in the commutator identity 



s=l 



(B.18) 



and apply the induction hypothesis. Finally, we trade i for p = i + s and use the identity 

1 



£ 



1 ^ "V-^,)*-^ 



f \ s ) \P ~ S J 1 — Q 



P. 



k- p (qd, 



'" l-q k 



(B.19) 



to simplify the result, proving (B.16) 



The next step is to rewrite (B.16) using bracket modes and Weierstrass functions. The 
key relation is 



EE 

i=0 fc^O 



hn-l + k\ 1 



t 



1-q 



1 x a t -h a +\ 



J2v m+1 (x,q)a[m], (B.20) 



m=0 



which is straightforward to derive by expanding and rearranging the sums. Differentiating 
with respect to r and integrating with respect to x, this also implies 



£ £ ( ha - 1 + k ) jfe-«tf / ! 



£=0 fc^O 



1-g* 



- 1 x a t -h a +i 



J29 i m+1 (x,q)a[m], (B.21) 



m=0 



where g l m+1 (x, q) is defined in appendix |A[ Setting a; = Zj/zi and summing over k in (B.16) 



completes the proof of the recursion relation (2.27) 
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C Details of the 6-point calculation 

In this appendix we provide additional technical details on the a 6 correction to the partition 
function outlined in section 13.31 Define 

Ilmk " ra = / it ' ' '/ if p ' +1 (I) Pm+1 W Pi+1 W Pr+1 (S) Pj+1 W ' 



Applying the recursion relation (2.27) repeatedly to (3.26) gives 



F((W, Zl ),...,(W,z 6 );r)= (C.2) 



2364J 
66561 



L2354J 
T66551 



75((W{£}W){r]d{k}W{m]W^)I £mrk [^''' 

+ 19(rf[r]^[A;]^[m]^[£]^_3)/, mfcr L2:!] , 

+ l20{(W[e]W)[r]W[0}W[k]W[m}W^)I £mkr [|1|] 

+ 14(W[r]d[ife]W[m]W[^]lV , _3>// rnrfc r66661 

+ 40((W[m]W[£}W)[r}d[k]W^)I emrk [ 23(<l j 

+ 30((iy[£]iy)[r]iy[A;]dH^-3)^r m [itli] 

+ 9(W[r]W[k]d[m]W[£]W- 3 )I ikrm [HI] 
+ 4(^[r]l^[A;]W[m]rf[£]^_3)/ mfcr , H 

+ 45((^^W^)[;]^[r]W[A;]^[m]W_3>i/r n fcri [fill 

+ 40((^[A;]^[m]^)[j]^[r]^[£]^_3)/£ mfcrj [f|l|] 
+ 5(W[;]W[r]W[A;]W[m]W[^W_ 3 )/ /TOfcrj - [S 
+ 15((iy[A;]W)[r](^HW / )[0]^[£]^_3)^ mfcr [f|ft] 
+ 30((W[£]W)L7'](W[m]W)[r]WMW_3>I^ri 



'634441 
.12536J 



where d = —W[0]W. We have discarded terms not linear in r because all such terms must 



cancel by the argument in section [3T2j The vevs are computed by applying the Jacobi identity 
(3.18) to remove nested bracket modes. Only the terms with bracket mode indices adding 
up to 5 will integrate to be finite, contour-independent contributions, so all others cancel. 
There are 37 non-zero contractions of the form 

(W[i]W[j]W[A;]W[/]W[TO]W-3} , (C.3) 

satisfying i+j + k + l + m = 5. These are computed with the help of Mathematica by 
applying the commutation relations and expanding out composite currents in products of 
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modes, giving for example, 

and other similar formulae. Each of the corresponding integrals can be computed straight- 
forwardly using the methods described in appendix |A| 



D The Woo [A] algebra 

Following [13] (see also [35]), but including some higher orders, the first few commutation 
relations of the the nonlinear W^ [A] algebra are 

[L m , L n ] = (m - n)L m+n + —m{m 2 - l)S m _„ (D.l) 



[L m ,W n ] = (2m-n)W n 



^3, ,, n , , , 



[W m , W n ] = 2(m-n)U m+n + — (m-n)(2m + 2n -mn-8)L m+n 



N 3 C / 2 iw 2 -u i 8 ^ 3 ^ uW 

——m[m -l)(m - 4)d m _ n H (m - njA^. 



[W m , C/„] = (3m - 2n)X m+n + -^(n 3 - 5m 3 - 3mn 2 + hm 2 n - 9n + l7m)W m+n 

i5iV 3 



25cN 3 [6m Zn > A ™+" 25cN 3 Um+n 



[U m} U n ] = 3(m - n)Y m+n - (m - n)^-^A)^ +n + n 44 (m - n)(m 2 - mn + n 2 - 7)U m+n 

CiV 3 

- —^-(108 - 39m 2 + 3m 4 + 20mn - 2m 3 n - 39n 2 + Am 2 n 2 - 2mn 3 + 3n A ) 
360 

cN 
x (m - n)L m+n - J^ m ( m2 ~ 1)(™ 2 - 4)(m 2 - 9)5 m _„ H 

[W m ,X n ] = (4m - 2n)y m+n - — ^(28m 3 - 21m 2 n + 9mn 2 - 2n 3 - 88m + 32n)U m+n 

56 N^ 

[U m , X n ] = (Am — 3n)Z m+n + n45(14m 3 — 5n 3 — 2lm 2 n + 15mn 2 — 86m + 59n)X m+n 

+ 5 — (28m 5 - 3n 5 - 35m 4 n + 30m 3 n 2 - 20m V + 10mn 4 

- 340m 3 + 355m 2 n - 220mn 2 + 75n 3 + 792m - A32n)W m+n + ■■■ 
[X m , X n ] = ^^m{m 2 - 1) (m 2 - 4) (m 2 - 9) (m 2 - 16)5 m ,_ n + • • • , 
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where L, W, U, X, Y, Z are the currents of spin 2, 3, 4, 5, 6, 7 respectively, and the constants 
are 

iV 3 = ^V(A 2 -4) iV 4 = -^a 4 (A 2 -4)(A 2 -9) 

5 6b 

N * = T^ * 6 ^ ~ 4 ^ " 9 ) ^ " 16 ) > ™ 44 = T-y^ - 19 ) 

105 15 

n 45 = J^°" 2 ( 3A2 " 97 ) n i = ^ 2 ( 29A2 - 284 ) • 

Here a is an arbitrary choice of normalization. We use the normalization of [2E], in which 
the WW OPE is 

W(z)W(0)^ 1 ^l + ... . (D.2) 

This corresponds to 

^ilA^ir < D - 3 > 

(For comparison to [26J, note that the charges are related by Wj ere = 27rWQ here ; for currents, 
however, jy here = 27T2W there , because the zero-mode in [26] is defined on the cylinder whereas 
ours is defined on the plane.) 



The dots in (D.l) indicate terms that will not contribute to the 6-point function of Ws 



at order c. Composite currents are defined as 

A (4) ~ LL , A (5) ~ LW , A (6) ~ WW , 6 (6) ~ -L'W + \LW' , (D.4) 

with the mode expansions 

A? = £ ■ L n- P L P : (D.5) 

A® = J2 ■ L n- P W P : 

A*? = Y, ■ W n- v W v : 

i 6) = £(|p - ») ■ L n- P W P : • 

These composite currents are not quasiprimary, but they can be made so by adding linear 
corrections. Since these linear corrections involve fewer fields, their contribution to the vevs 
will be subleading for large c, and hence do not affect our results. 
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